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1 Introduction 



Two-dimensional quantum gravity is important to explore the nonperturbative dy- 
namics of string models and statistical models of random surface as well as to provide 
a theoretical laboratory for quantum gravity in higher dimensions. The conformal 
field theory approach has successfully described two-dimensional quantum gravity 
and evaluated its universal quantities [Q. However, the continuum description has 
encountered a difficulty where quantum gravity is coupled to the conformal matter 
field with the central charge c > 1, despite its importance furnishing off-critical 
string models in an arbitrary dimensional target space. 

As a discretized approach, the matrix model 0, |3| provides a more rigorous 
and powerful formulation of two-dimensional quantum gravity and off-critical string 
models. Exact solutions in matrix models have been obtained so far only for cases 
corresponding to matter theories with c < 1 coupled to quantum gravity It 
is straightforward to define matrix model candidates for quantum gravity coupled 
to conformal matter with c > 1, such as multi-Ising models 0. However, it has 
been difficult to solve these theories exactly. Therefore it is most desirable to find 
a reliable approximation method which gives correct results for solved models and 
reasonable results for so far unsolved models. 

Brezin and Zinn- Justin have proposed a large- iV renormalization group (RG) 
equation as such an approximation method || Subsequently, we have explicitly 
derived exact RG equations for 0(iV)-vector models |T0| , llfl and one-matrix models 



I2| . We found that it is crucial to take account of the reparametrization identities 
in order to obtain a meaningful RG equation. Moreover, the RG equation turns out 
to be a nonlinear equation for the matrix model in contrast to the linear one for 
vector models. 

The purpose of the present paper is to explain our large- N RG approach for 
matrix models in depth and to report on results for a two-matrix model. We obtain 
correct locations of several nontrivial fixed points and scaling exponents. With 
this result, we can identify the matrix model representation of the gravitationally 
dressed operators corresponding to those in the continuum theory. Moreover, we 
find that the large- N RG is useful to describe the global phase structure of matrix 
models. The global picture of the RG flow can be drawn practically by a linearized 
approximation for the nonlinear RG equation. It should be noted that the linearized 
approximation is meaningful only after the reparametrization identities are taken 
into account. We observe from this global flow that an RG trajectory connects a 
pair of fixed points corresponding to two unitary conformal matter theories coupled 
to gravity, in accordance with the expectation based on conformal field theories 
over fixed background This flow can be understood by means of the usual 

Kadomtsev-Petviashvili (KP) flow [fP|]-||16|. On the other hand, we also find that 
a pair of fixed points are connected by an RG trajectory which is out of the KP 
hierarchical description. From these examples, we conjecture that a kind of the 
gravitational analogue of Zamolodchikov's c-theorem holds in accordance with some 
previous expectations |T7| for several cases. 



The partition function Z^(gj) of a one-matrix model with a general potential 



V{<fi) is defined by an integral over aniVxiV hermitian matrix <p 



Z N { 9j ) = j #exp[-jVtrV(0)], (1.1) 
k>i K 

In particular, we are interested in the subspace {g± — 0, g-i — 1, gj = for j > 3 m}. 
Let us define the free energy by 



Z n(9j) 



Zn(92 = 1, others = 0) 



;i-2) 



The free energy of the matrix model gives a generating function of connected Feyn- 
man diagrams constructed in terms of j-point vertices with the coupling gj. By con- 
sidering the dual of the Feynman diagrams, we can interpret that the dual diagram 
gives a decomposition of a two-dimensional surface in terms of j-gons corresponding 
to the j-point vertices gj. Under this interpretation, the free energy of the ma- 
trix model is identified with the partition function of the two-dimensional quantum 
gravity. One finds that contributions from the different genera are distinguished by 
powers of 1/N 2 . 

For brevity, we consider the triangulation (g 3 = g, gj = for j > 3). The 
free energy F h (g) for genus h surfaces is given by the coefficient of N~ 2h in the 
1/iV-expansion, and the number n of triangles is counted by the power of g 

oo 

F(N,g) = J2N- 2h F h (g), F h (g) = 5>V fc (n). (1.3) 

h=0 n 

Here Vh{ n ) is the number of Feynman diagrams with n vertices and genus h divided 
by the symmetry factors. The free energy F h (g) can be regarded as the partition 
function of the randomly triangulated surface with genus h. Its singular part behaves 
as 

-^skig(50 ~ (9 ~ 9*) ~ <T,i_2 ' ) + less singular terms (1.4) 

around a critical point g*. The singular behavior (|1.4j ) implies that the number of 
triangles Uh(n) has an asymptotic form 

v h {n) ~ n Th - 3 g; n (n -> 00). (1.5) 

The average number of triangles (n) of surfaces with genus h is estimated by ( |1.4| ) 

n) = — — 

En 9 M n > 

= g |-logFi g ( g )~ (r ;~ 2 ^ ^oo (g^g m ). (1.6) 
^9 9* 9 

The average number of triangles in triangulations with any genus diverges simulta- 
neously at g — g* with the same exponent —1. 

In matrix models (and vector models fl8|) solved so far, the susceptibility expo- 
nent depends linearly upon h as 

I\ = 70 + hji. (1.7) 



Therefore we can take the double scaling limit 



l/N - 0, 



9 -> 



with N~ 2hl (g - g. 



fixed, 



and retain nontrivial contributions from all genera to the free energy. In this limit 
the singular part of the free energy takes the form 



F sing (N,g) = (g- g ^f(N 2 ^{g-g*)), 



1.9) 



where / is an unspecified function. 

Let us discuss the double scaling limit ( |1.8|) from the viewpoint of the continuum 
limit in two-dimensional discretized quantum gravity. We consider the following 
continuum limit for each genus h by sending the length a of edges of triangles 
(spacing of the random lattice) to zero. If the average number of triangles (n) 
diverges, we can scale the length to zero so as to keep the 'physical' area fixed: 



0. 



<n> 



oo, 



with 



(n) = fixed. 



;i.io) 



Since the scaling behavior ( |1.6| ) shows that the average number of triangles diverges 
by letting g approach a critical value g*, we can regard the double scaling limit ( |1.8|) 
as the continuum limit eq.( |l.ld| ) for each genus h 

N~ 2hl ■ (n) = fixed a 2 ■ (n) = fixed. (1.11) 

Therefore the double scaling limit suggests that jV _1//71 is proportional to the lattice 



spacing a. A similar point of view has been expressed previously in refs. M, IT9 



Brezin and Zinn- Justin proposed the following approximation scheme for calcu- 
lating the location and the exponents of critical points [|]. Integrating out parts of 
(N + 1) x (TV + 1) components of the matrix field in eq. (|1.1| ) gives us the following 
recursion relation in the first nontrivial order of perturbation: 

Z N+l {g) = \(N,g) N2 Z N (g + 5g(N,g)). (1.12) 

One expects that eq.( |1.12| ) is well approximated by a linear differential equation 



N 



_d_ 

dN 



P(9)§- g +l(9) 



F(N,g) = r(g). 



;i.i3) 



The homogeneous part (i.e. r(g) set to zero) of eq.( |1.13D determines the singular 
part of F around a zero g = g* (a fixed point) of the /^-function as 



F smg (N,g) = (g-g^*)m 9 *)f(N^*\g-g*)). 



;i.i4) 



This singular behavior is just the same as that of the exact solutions (|1.9| ) of the dou- 
ble scaled matrix model. Namely, the string susceptibility exponents are identified 

as 



7o + fryi 



;i.i5) 



Therefore once one knows the functions [3(g) and j{g), one can compute positions 
and exponents of critical points. In fact, the RG equation computed in the one- 
loop approximation provides us with a fairly good result || [2(J. However we have 
confirmed in ref . |T]J that no improvement is achieved if one pursues the perturbation 



to higher-loops by simply truncating induced interactions of higher orders. We now 
understand that the naive perturbative method fails due to the reparametrization 
invariance of the partition function, as explained in sect. 2. 

Eq. (|1.13 ) resembles RG equations in field theories in its form. Since eq. Ql.llj) 



shows that we can identify N^ 1 ^ 1 with the lattice spacing a, we can regard an 
infinitesimal shift of N as a substitute for a shift of the lattice spacing a. Thus 
it may be appropriate to call a differential equation describing the response of the 
physical quantity under the infinitesimal change of N, such as eq. Q1.13Q , an RG 
equation for matrix models. 

This paper is organized as follows. In sect. 2, we derive exact RG equations 
for one- and two-matrix models via coset and/or eigenvalue methods. In sect. 3, 
we develop a scheme to solve the RG equations and show that they reproduce the 
spectrum of the exact solutions. In sect. 4, we analyze the topology of the RG flow 
using the linear approximation to the RG equations, and discuss the gravitational 
c-theorem. Possible applications of our method, including generalizations to can- 
didates for c > 1 quantum gravity, are addressed in sect. 5. The reparametrization 
identities for the two-matrix model are solved in appendix A. 



2 Large- N renormalization group equation 



2.1 One- matrix model I: coset method 

RG transformation 

In this subsection, we restrict ourselves to the simplest case m = 3 and set 
g = g 3 > without loss of generality. 

To derive a relation (|1.12|) , we employ a decomposition of an (N + 1) x (N + 1) 
hermitian matrix variable $ into an N x N hermitian matrix <f), a complex iV-vector 
v, its transposed conjugate and a real scalar a: 






V 




a t 



(2.1) 



In terms of <ft, v, v* and a, the partition function ( |1 . 1|) reads 
Z N+l (g)= I d(j)dvdv ] da exp \—(N + l)(trV(0) + V(a) + v\l + g((f) + al))v) 



First we integrate over v and exactly to obtain 

Z n+M = (VrVT /#exp[-(iV + l)trV(0)] 



(2.2) 



N + l 



J daexp[-(N + l)V(a) - tr log(l + g(<j> + al))). (2.3) 



Now we are ready to evaluate the a-integral by the saddle point method, system- 
atically in each order of the 1/iV-expansion. The saddle point a s = a s (g,(f)) is 
determined by the saddle point equation 



2 , 9 ... 1 



°- + g<+AP tr l + ^ + a ,l) =0 (2 ' 4) 
as a U(N) -invariant function in <fr (i.e. a function in tr0 J, s). 

If we explicitly retain only the leading part of the 1/iV-expansion, we have 
Z N+ i(g) ( vr ^ 



z N{9) kn+u M-tovw-Nvw 

- tr log(l + g{cj> + a s l)) + 0(N°)}\ . (2.5) 



Here (■ • •) denotes the average Z^ 1 J d(j){- • -) exp[— iVtr V(^)]. In the following we 
repeatedly utilize the factorization property of a multi-point function of U(N)- 
invariants O, O' into a product of one-point functions, 

(O O') = (O) (O 1 ) + 0(N~ 2 ) (2.6) 



in the large- N limit. In view of this factorization property, eq.( |2.5|) can further be 
written as 

z N+ i(g) ( ti \ N 
Z N {g) \N+h 

- ( tr log(l + g(<f> + (a s ) 1))) + 0(N°)] . (2.7) 



exp[-(trF(0))- NV{{a s )) 



Finally we approximate a difference in eq.( |2.7| ) by a differential to obtain 

F(N,g) (2.8) 



ON 



= ~\ + tr + V((a s )) + ^tx log (1 + g{<t> + (a s ) 1))) + O (1 

Reparametrization invariance 



The right hand side of eq.( |2.8| ) is a function of variables N, g and ((1/N) tr0 J ) = 
j(dF/dgj) (j = 1, 2, . . .). What we will do here is to express the dF/dgj- dependence 
(J 3) in terms of dF/dg^. This may be done by taking into account reparametriza- 
tion invariance of the partition function (|1 . 1|) . In this way, we can reduce all the 
induced interactions to those in the original potential. 



For a while, we work on a generic potential (|1 . 1|) . The partition function (|1.1|) is 
expressed as an integral and we can freely reparametrize the integration variable </> 
without changing the partition function. Obviously only reparametrizations which 
are analytic at the origin in the field space are allowed. If we perform reparametriza- 
tions 4>' = (ft + e(j) n+1 (n > —1) for the coupling constants receive shifts. This 
reparametrization invariance implies identities 

" \N 



i=0 



m , 1 

E^(T7tr0" +fc ) (n>-l). (2.9) 



k=l 



We note that eq.( |2.9| ) constitutes the first (K = 0) set of full reparametrization 
identities (discrete Schwinger-Dyson equations) Jl6], |21| 

/#tr-f- [<p n+l TT tr0 nfc exp[-iVtr1/(0)]) = 0, (n > -1). (2.10) 
J d<p { ^ 1 J 

Only reparametrization identities in the large- iV limit will turn out to be necessary 
in our context. When we truncate the full set of identities (|2.10| ) to those in the 
large- N limit, eq.( [2.9|) exhausts all independent relations, due to the factorization 
property ( |2.6| ). 



Eq.( |2.9| ) can be neatly expressed by introducing the resolvent operator g] 

The expectation value (W(z)) is the generating function of one-point functions: 

00 /I \ 

(W(z)>=5:z- 1 -V-tr0»). (2.12) 

3=0 XiV 1 

Using the resolvent, eq.(|2.9|) is summarized into a single equality for \W(z] 



(W(z)) 2 - V'(z) (W(z)) + Q(z) = 0, (2.13) 

-i^g<^(V)G*'>- (-» 

Eq. fl2.13|) is usually referred to as the loop equation. In the above, Q(z) consists of 
(ticf) k ) for 1 < k < m — 2. For our purpose, one-point functions (tr0 fc ) for 3 < k < m 
are to be regarded as independent variables, corresponding to the coupling constants 
present in the original potential. Therefore a special treatment must be applied for 
(tr<j)) and (tr</> 2 ); we should replace them in terms of (ti(j) k ) for 3 < k < m with the 
help of the first two equations of 



= (Itr^ + f^/ltr^ 1 ), (2.15a) 



\N V ts \N 
1 = (^tr0 2 ) + f:^ fc (ltr0 fc ). (2.15b) 

Now let us denote the one-point function of a single trace as aj = ((1/A^)tr0- J ), 
and that of the resolvent with the above replacement as W. It satisfies 

W(z; g 3 - a,) 2 - V(z)W(z; 9j ; a,) + Q(z; 9j - a 3 ) = 0, (2.16) 

with Q(z; gj; aj) given by 



Q( z i 93: ■ ■ ■ ? 9m] • • • j o TO ) 



m— 1 



= E^ (fc+1) W(- 
fe=i 

m— 1 

fc=2 
m— 1 

+ E^ (fc+1) W(- 



\fc-i 



jfe! 



t, _ i m m 

V=4 V £=3 



\fc-l-2 



fc! 

1 (Jfe -!)(*- 2) 



fc=3 

m— 2 m— 1 



+ E E v^\z){- 

j=3 k=j+l 



2 ■ 2 



i-E 

e=3 



k\ 



jfe - 1 
J 



j a 



j- 



(2.17) 



The quadratic nature of the loop equation (|2.16 ) is a peculiarity of the one-matrix 
models, as compared to the multi- matrix models. For the case m — 3, the explicit 
form of W(z, g,a = a^) reads 



W{z,g,a) 



2\ Z + 9 Z 



> + ^ 2 ) 2 - 4(1 + #2 - # 2 + 3# 3 a)) . 



(2.18) 



Note that, by definition Q2.11| ), only a solution with the large- ,2 asymptotic behavior 
W (z) ~ 1/z is allowed. Moreover, the eigenvalue density distribution p(z), related 
to W(z) by p(z) — — l/V Imiy(z), should satisfy the trivial identity 



p(z)dz = 1. 



(2.19) 



We have observed that these two requirements are enough to determine without 
ambiguity for the present well as in the case of the two-matrix model. 

Another simplification occurs for the averaged form of the saddle point equation 



(a s ) + g {a s )' 



W I— - (a B ),g,a 



0. 



(2.20) 



By combining it with the loop equation ( ^.1G| ), we explicitly obtain (a s ) in the form 



-g + 3g 2 a = a(g,a). 



(2.21) 



This is in fact an expected result, because (a s ) = (& S N+1N+1 ) should be equal to 

(jjtT<p\ = (^(1/N) J2iLi This fact helps computation in the coset decomposi- 
tion method considerably, for the case of two-matrix model as well. 

Using eqs.( p78|) , Q2.18Q and (|2.21| ), we obtain a differential equation obeyed by 
the free energy F(N,g) of the one-matrix model with cubic coupling, 



d 

N 0N + 2 



F(N,g)=G(g,%-\ +0(-). 



dg 



N 



(2.22) 



G(g, a) 



-a+-a(g,a) +-a(g,a) + log (1 + ga(g, a)) 



+ 



dz ( W(z, g, a) 



Since the above equation describes the response under the shift iV — > A" + 1 it de- 
serves the name of RG equation, although it is (non-polynomially) nonlinear unlike 
the form anticipated in eq. (|1.13|) . 



Comment on the nonlinearity 

One might wonder if we can keep the equation linear at the cost of having 
an infinite dimensional coupling constant space. We could generalize the partition 
function to include interactions such as tr<p ni tr<ft n2 ■ ■ -, 

Z N {g ni n 2 -) = J ' d(j> exp \-NS{(j>)} , (2.23) 

oo k 

S(<t>) = E E g ni n 2 ...n k N l ~ k Y[n7Hr^- 

k=l (m,n2,...,n k ) i=l 

and enlarge the coupling constant space to {(gi, <?2, • • • , gn,gi2, ■ ■ ■ , gn^n 2 --n k -, ■■■)}■ 
Then we would expand the right hand side of eq. (|2.8| ) in terms of 




(N, 0,1, g, 0,0,...) (2.24) 



to obtain a linear 'RG equation' of the form 

dF 

F(N, 0, 1, g,0,0,...)= E Pm-nM* 0, 1, g, 0, 0, . . .). 

(m,...,n k ) °9m-n k 

(2.25) 

It is straightforward to check that there exists no simultaneous zero for each function 
f3 ni ... nk (g). This contradicts the exact solution of the one-matrix model, which has 
a nontrivial fixed point at g\ — 0, g^ = 1,<73 = 432 _1//4 , gj = (j > 4). The obvious 
reason for this apparent contradiction is that the (dF/9g ni ... nfe )|( 9j ) = (o,i i9i o,o,...) are n °t 
independent quantities in the large- N limit. Rather, they are mutually dependent 
through an infinite number of relations ( p.6|) and ( |2.9| ) and are therefore '/3-functions' 
for g ni ... nk . Accordingly we are forced to take the reparametrization identities into 
account to exclude this ambiguity of /5-functions. 



d 

N m + 2 



2.2 One-matrix model II: eigenvalue method 

The RG transformation by means of the coset decomposition is in principle applica- 
ble for an arbitrary potential as well. We can eliminate terms such as (v >v ) 2 , which 
are present in the quartic potential case, by introducing auxiliary fields. However, for 
the model which admits an eigenvalue representation via U (AQ-integration we have 
a more convenient RG transformation: integration over the (N + l)-th eigenvalue, 
the extra degree of freedom. This procedure is particularly useful for one-matrix 
models with a generic potential (m > 4), since it can be applied to these models on 
the same footing. 

We can integrate over the U(N) variables in ( |1.1|) to obtain an integral over the 



eigenvalues {Xj} [0 

N 

ZN(gj) = c N J Y[d\i A N (\i) 2 exp 



i=l 



iV 



i=l 



(2.26) 



where A N denotes the Van der Monde determinant A N (\j) = rii<i<j<w(^i ~ -\?) 
and the volume of the U(N) group = ■K N ( N ~ l )i 2 j J\^ =l p\. 

In order to relate Z^+i with Zn, we shall integrate over the (iV + l)-th eigenvalue 
A = Atv + i in Zjy + i 

Z N+1 ( 9j ) = jd® exp[-(JV+l)try($)] (2.27) 

. N w . N 

= c N+1 [ Y[dh AliX,) e-^ +1 )E i=1 ^)| dAn | A _ Ai |2 e -(iv + i)y(A) 

i=l ' i=l 

I ' d<p e -(N+i)*v&) J dX eX p[-(A^+ l)F(A) +2tr log |A1 - 0|] . 



CN+l 



C N 



The A-integral can be evaluated by the saddle point method as a power series in 1/N 
around the saddle point, since the effective potential (N + 1) V(A) — 2 tr log | Al — <j)\ 
is of order 0(N l ). The saddle point equation 



V'{\ s ) 



2 1 
iV tr AJ^ 



OO 1 

,-n-l 1 



2£A S 

n=0 



N 



tr, 



(2.28) 



implicitly determines the saddle point A s = X s (gj,<p) as a ?7(iV)-invariant function 
in 0. By substituting A s into the partition function, we find 



2jv(5j) 



(exp [- tr V(0) - iVF (A 9 ) + 2 tr log I A s l -<f>\ + 0(N°)] ) (2.29) 



exp 



( tr V{(j>)) - NV{{X a )) + 2 ( tr log | (A s ) 1 - 0|> + 0{N 



0> 



Here the factorization property (p.6|) is used in the second line, and (A s ) = (X s (9j, (ft)) 



is determined again using the resolvent p. 11 



n<A s » = 2(ltr. ' 



N (A s )l- 



(2.30) 



By taking the large-iV limit, we get the following differential equation as an RG 
equation for the free energy, 



F(N, 9j ) (2.31) 
+ tr V {<!>)) + V((X S )) - 2 tr log |(A S ) 1 - 0|) + O (^) . 



As in the case of the coset method, we can express the right hand side of eq. ( 2.31]) 
in terms of gj and aj = dF/dgj for j — 3, • • • , m with the help of reparametrization 
identities. It is straightforward to express (A s ) in terms of these quantities, since the 
saddle point equation ( |2.30|) readily takes the form V'((X S )) = 2W((X S ) ; gj; a,j). By 



combining it with the loop equation ( 2.16| ), we find an equivalent and more useful 
expression 

H(As» 2 -4Q((A s >;^; % ) = 0. (2.32) 
Eq. fl2.32D shows that (A s ) falls on one of the edges of the eigenvalue distribution 



p(A), because the left hand side is exactly the descriminant of the loop equation 
(|2.16|) . This fact can be understood as follows. We adopt the picture that the 
N + 1 eigenvalues, confined by the potential V^(A), interact with each other via 
a repulsive Coulomb potential log|Aj — Xj\. Consider the effective potential for 
the (N + l)-th eigenvalue, which is generated by other N eigenvalues obeying the 
distribution function p(A). If the (N + l)-th eigenvalue falls on a point p(A) ^ 0, 
it costs an energy loss (of order N~ l log N compared to the total energy) from the 
Coulomb interaction. On the other hand, if the total action is to be minimized 
when the (N + l)-th eigenvalue is placed on an isolated point outside the support of 
p(A), the original distribution is unlikely to minimize the action for N eigenvalues. 
Consequently the (N + l)-th eigenvalue should fall just on one of the edges of 
p(A). Moreover, the true saddle point should converge in the g — > limit to 2 
(Wigner distribution), and among such roots it should minimize the action S(X) = 
V(X) — ((2/iV)tr log |A1 — We can always select the true saddle point out of 
the 2m — 2 roots of eq. (|2.32|) by imposing these two requirements, and denote it as 
A = X(gf, ai). 

In this way, we can express the right hand side of eq. ( |2.31| ) in terms of gi and 
dF/dgi (i = 3, . . . ,m), 



n Wm + 2 



CdF dF \ ( 1 



G r eig(5 , 3, ■ ■ ■ ) 9m] ^3; • • • j a m) (2.33) 
m k — 2 

= -l-J2^-9ka k + V(X)-2logX 

k=3 1 



-2 [ dz (w(z; gf 3 , . . . , g m ] a 3 , . . . , a m ) - - 

J±oo \ Z 



The complete set of nonlinear RG equations for the one-matrix model consists of the 
above equations fl2.33j ), ( [2.16| ), and (|2.32|) . So far we have much numerical evidence 



for the equality of the function G in the coset method and the function G e i g in the 
eigenvalue method, but have not been successful in showing this analytically. 

2.3 Two-matrix model 

We consider a two-matrix model defined by 

Z N (g + ,g_,c) = J d0 + d0_exp[-iVtr(y(0 + ,0 + ) + y(0_,p_) + c0 + 0_)], 

V{x,g) = ^x 2 + 9 -x\ (2.34) 

where <p± are N x N hermitian matrices. The free energy 

771 /at \ 1 i Z N (g + ,g^,c) rt „ rt 

F(N, g + , c) = -- log Zn{qaq) (2.35) 



represents the partition function of the Ising model on random triangulated surfaces 
23] , and the Ising temperature T is proportional to — log _1 (— c). This model is also 



known to allow an eigenvalue representation [24j 

2\jV(AT-l)/2 „ n 



Zn(9+,9-,c) 



-CTC 



nfp! 



i=l 



exp 



N 



-N^{V(X+ dig+ ) + V(\-j,g-) + cX+jX-j) 

3=1 



. (2.36) 



However, the integrand is an alternating function in A-^j, since the Van der Monde 
determinants is not squared but bi-linear in the above expression. Therefore it is 
inappropriate to evaluate the A± 5 Ar + i-integral by the saddle point method. This argu- 
ment also applies to multi- matrix models with a chain-like interaction c Y?a=\ ^ T (fta^a+i- 
For this reason we present only the coset method in the following. 



RG transformation 

We decompose (N + 1) x (TV + 1) hermitian matrices $± into 



' 4>± 


v± 


, 4 


a ±/ 



(2.37) 



as with the case of one-matrix models. In terms of these variables, the partition 
function (|2.34j) reads 



Zn+i{9+,9-,c) 



4> + d(p-dv + dv\_dv-dv^_da + da- 

■ exp[-(iV + 1) tr (V(<j>+, g+) + V(4>-,gJ) + c0 + 0_)] 
• exp[-(JV + l)(V(a + , g+) + V{a.,gJ) + c«+a_)] 



(2.38) 



• exp 



■(N + l) ■ (vlvl) 



1 + ^ 



cl 



cl 



1 + 9- 



a-1) 



The part of the action quadratic in (v + , u_) and (v\_, vi.) is presented in the form of 
a 2N x 2N matrix. We perform the integration over (t> + , t>_)-field to obtain 



ZN+i(g+,g-,c) 



7T 



N+l 



2iV 



(2.39) 



■ exp [-(N + 1) tr (V(<f>+, g+) + V(<f>.,g-) + c0 + 0_)] 

• J da + da_ exp[— (N + l)(V(a + , g + ) + V(ct-, gJ) + ca + a_)] 

1 + 9+(4>+ +a+l) cl 

cl 1 +g-U- + a_l) 



■ exp 



-Tr log 



The trace Tr is taken over a 2iV x 2iV matrix, whereas tr is over an N x N matrix. 
Again we employ the saddle point method for the a±-integration. The saddle point 



a S) ± = a St ±(g+, g~, </>+, </>_) is determined by the equations 



-i 



«s, ± +c« SiT +^a s , ± +^-Tr ± ^ cl i + ^(^_ + Qbi _1) ) = °' 

(2.40) 

where 

Tr ± M = Tr (P±M), P + = ( J °V P_ = ( ° jV (2.41) 

As in the case of the one-matrix model, we employ the factorization property 
and approximate Z N+1 /Z N by a differential to obtain 



F(N,g + ,g.,c) (2.42) 



= -1 + (^ tr ^+> 9+)) + (^teV(<f>-,9-)) + c ^tr<M>-) 
+V'«a Bi+ ) , 5-+) + V^Og ,_) , flf_) + c (a Si+ ) (a 8 _) 

+ /l T r log 1 + ^ + + 1} Cl "| \ + O (- 

where 

(•■■) = ^ / ■ ■) exp[-iVtr (V(</> + , g + ) + 1/(0-, <?-) + c0 + 0_)]. (2.43) 

Reparametrization invariance 

In what follows we concentrate on the case g + = g_ = g which describes the Ising 
model coupled to two-dimensional quantum gravity in the absence of an external 
magnetic field. In addition, we assume that 

(ti </>*+) = (tr <f>i) (2.44) 

holds for any i. When we trace on a part (T < T cr ) of the critical line, this as- 
sumption amounts to defining the free energy as the sum of contributions from two 
spontaneously magnetized states. A direct consequence of this assumption is 

(a Si+ ) = (a s ,_> = (a s ) . (2.45) 

Under these conditions, the saddle point equation for (a s ) reduces to the form 

(« s ) + c (a.) + g (« s ) 2 + /^Tr 1 — - \ = 0, (2.46) 

\2N (l/g+ (a s ))cr + (c/g)tri + 0/ 

where 

(i!)-"- (!!)■ (2A7) 

The differential equation ( 2.42[ ) simplifies into 



F(N,g,c) 



= -l + (^trV(0 + ,(7)^ + ^trr(0_,^)^ + c^tr0 + 0_^ (2.48) 
+ 2V((a s ) ,g) + c (a s ) 2 + /^Tr log ( (- + (a s )j <r + -a x + <f> J \ + O (-^ 



The right hand side of eq.( 2.4£p consists of terms of the form tr (0" 1 0™ 1 0" 2 0!! 12 • • •) 
We observe that all these induced interactions can be reduced solely to the interac- 
tions present in the original potential, i. e. tr (0+ + and tr by recursively 
using the reparametrization identities 



tr- 



d 



{0+ 1 0™ 1 0+ 2 0_ 2 • • • 

• exp[-iVtr (V{<f>+,g) + V(<f>-,g)+af> + <l>-)]} = 0. (2.49) 

Again this procedure can be performed at once by introducing the resolvent operator; 
in fact we can find a closed set of loop equations (A. 2) that expresses the one-point 
function of the resolvent of the required form (see the appendix A) 



W (z;g,c; a g ,a c 



:Tr 



N za + (c/g)<r 1 + (j> 



in terms of a g = ((1/3JV) tr (0% + and a c = ((l/N) tr <p+(p-). 
Furthermore, by combining the saddle point equation (|2.46| ), or 



(2.50) 



1 



1 



(1 + c) (a s ) + g (a s ) + -W ^- + (a B ) ; g, c; a g , a c 
and the loop equation ( |A.4| ) we can prove 



(as) = ( y ' »' V 



9 



1 + c 



■l+ca c + ) = <>(.'/• r: a ir ". ' 



(2.51) 



(2.52) 



as expected. These three equations (|2.48|) , ( |A.4| ) and (|2.52|) constitute the nonlinear 
RG equation for the two-matrix model 



N w + 2 



„ „ I dF dF \ / 1 



(2.53) 



G(g, c; a g , a c 



9, , /i , _x=-2 , 2 _ 3 



|a 9 + (1 + c)a 2 + -3a 3 + 2 log(l + <?a) 
+ / dz [W (z;g,c;a g ,a c ) . 

J— oo V ZJ 



3 Fixed points and operator contents 
3.1 Solution to the nonlinear RG equation 

In the previous section we have derived nonlinear differential equations, governing 
the critical behavior of the free energy F(N, g), of the form 



d 

N 0N +2 



F{N,g) = G[g, 



dF(N,gY 
dg 



(3.1) 



(or its multi- variable extension) instead of the linear one ( |1.13[ ). To complete the 
program we now need a formula which relates positions and exponents of critical 
points to G, just as in the linear case eq. (|1.15| ). 



Sphere limit 

First we concentrate on the leading part F°(g) of the free energy in the 1/N- 
expansion. It is easy to see that F° satisfies 

dF°(g)\ 



2F°(g) = G (g, 



dg J 



(3.2) 



We assume that F° consists of analytic and non-analytic parts around a critical 
point 

oo oo 

F (9) = T,^(9-9,) k + T,U9-9*) k+2 ~"° (7o * Z). (3.3) 

fc=0 fc=0 

Quantities a^, 7o, g* are unknown and are to be determined. In anticipation of the 
result, we have used the same notation as the susceptibility exponent 70 in eqs.( |1.4| ) 
and (|T?T|). We expand the function G(g,a) around (g*,ai) as follows 

00 

G(g,a) = Y,f 3 n(g)(a-a 1 ) n 

n=0 

00 00 

= E £ A* -to -oi) B - ( 3 - 4 ) 

n=0 fc=0 

We substitute these expressions into the nonlinear RG equation ( |3.2|) and compare 
the coefficients of various powers of g — g* on both sides. Since b 7^ by definition, 
the most singular term (g — g* ) 1-70 determines the fixed point 

= p 10 (2 - 7o)&o- (3.5a) 

The next singular term (g — g* ) 2-70 determines the critical exponent 70 

2b = (2 - 70) W11 + &o ■ 4o 2 ). (3.5b) 

So far there appear two unknown coefficients a\ and a 2 from the analytic part. 
However, by comparing terms (g — g*) and (g — g*) 2 , we can fix these quantities 

2ai = Poi, (3.5c) 
2a 2 = A)2 + Ai2a 2 + /3 20 (2a 2 ) 2 . (3.5d) 

These four equations can be rewritten in terms of the function G(g, a): 

= G a (#*,ai), (3.6a) 

2 

7, = G, g ,a(g*,a>i) + 2a 2 G>,a(#*,ai), (3.6b) 

* — 7o 

2a x = G )B (g*,ai), (3.6c) 

2a 2 = -Gg,^*, ai) + 2a 2 G iSi a(^, a x ) + 2(a 2 ) 2 G^ a {g^ a t ). (3.6d) 



Obviously eqs. ( |3.6a| ) and ( |3.6b[ ) are the nonlinear counterparts of the usual equations 
determining fixed points and scaling exponents, 



= /%>), and 2/(2- To ) = /3 , Gfc), (3.7) 

respectively. We can solve the above four equations simultaneously to fix four un- 
knowns, the critical coupling g*, critical exponent 70, and coefficients ai,a 2 . More- 
over, the other consistency conditions determine all the other coefficients at each 
fixed point, a k (k > 0) and b k /bo (k > 1) recursively except for the overall normal- 
ization of the singular term b . Although the equation for a 2 is quadratic, we can 
choose the correct branch, that which continues to the unique solution at the origin 
(the Gaussian fixed point). In this way we can obtain the series expansion of the 
sphere free energy around the fixed point, up to one integration constant. 

It is straightforward to generalize our algorithm to the case of multi-coupling 
constants. We consider a differential equation 



NjL + 2 
ON 



maHG ( ft; g + o(i). p. 



We assume that the leading part F°(N,gj) of the 1/iV-expansion of a solution can 
be expanded in the form (Sgj = g,j — g^) 



F°(N, g 5 ) = J2 a h-jJ9h ■ ■ ■ 5g jn + ( v ij S 9jf 70 b (i)h-jJ9h ■ ■ ■ Sg jn , (3.9) 



n=0 i n=0 



where the repeated indices are summed over. A transformation matrix is in- 
troduced to diagonalize the scaling matrix. Then we have a set of equations to 
determine g^, % \ dj and 

dG 

= _, (3.10a) 

T 2-7^ ] dgjdai daida k 

dC 

2a t = — (3.10c) 

% 

1 d 2 G d 2 G d 2 G n d 2 G 

2a v = o^ - ^ h a ^^~^ h a ik^-^ H 2a ik aji- — — , (3.10d) 

2 ogiOgj ogjOa k ogiOa k oa k oae 

where they stem from terms of order O (^(VijSgjY'^o ^ , O (^Yi k 5 gk) 1 ^ 1 ^ {VjkSgSj > 

0(8 gi) and 0(5gi5gj), respectively. We refer to the right hand side Gij of eq. (|3.10b| ) 
as the 'scaling exponent matrix'. 

Exceptional cases 

So far we have implicitly assumed that the exponent 70 is irrational. In most cases 
this assumption is incorrect and higher powers of -F S i ng appearing in the right hand 
side of eq. (|3~2j) contribute to the analytic part. Therefore eqs.( |3.6c| -d) which serve 



to determine the coefficients ai and a 2 may require modification on such occasions. 



However, for c < 1 theories that we examine in this paper, 70 is known to be 
negative. Then we observe that the potentially dangerous terms 

(dF sing /dg) k ~ 0(5g^ k ) (k > 2) (3.11) 

are always of higher order than 0(5 g 2 ). Hence the mixing of such terms does not 
occur in eqs. (|3.6cj -d) determining a± and ci2 which are needed to identify fixed points 
and scaling exponents. 

The function /3 n (g) in eq. (|3.4j) may sometimes develops singularity. In such oc- 
casions the G-function can be separated into its regular and singular parts around 
a specific point according to the power expansion of each /3 n (g). Then we can prove 
that our algorithm, which originally assumes the analyticity of the G-function, is still 
valid after a replacement of G in eqs.( p.6a| ~d) or ( |3.10a| -d) with its regular part C7 reg . 



Higher genus contributions 



Now we show that eq. ( |3.1| ) is sufficient to determine the singular behavior of 



higher genus free energies F . Suppose we retain all terms of subleading order 
in 1/N in each procedure, i. e. saddle point evaluation of the a- or A-integral, 
application of the reparametrization identities, and replacement of a difference in iV 
by differentials. Then we can express this complete RG equation in the form 



F(N,g) = G (g,?P\ + AG(N,g;[F\), (3.12) 



dg J 

where [F] denotes differential polynomials of F(N,g). 

The partial differential equation ( p.!2j ) can be separated into a set of ordinary 
differential equations for each genus contribution. It is important to realize that the 
additional contribution AG introduced into the right hand side of the RG equation 
carry additional powers of 1/N. Therefore they will not appear in the coefficient 

00 n rfphi 

{2-2h)F\g)=r\g) + Y.[3n(g) £ H -r-. (3.13) 

n=1 /!. 1 +--- + /ln=/l,*= 1 " 

hi>0 

The inhomogeneous term r h (g) does get various higher order contributions including 
those terms from F h for h' < h — 1. 

Now we assume that F h {h! < h — 1) is less singular than F h . Then the leading 
singular part in the right hand side of eq.( |3.13j ) is dominated by the terms involving 



dF h /dg. Once we realize this structure of the RG equation, we can repeat the same 
argument as in the case of the sphere to obtain the scaling behavior of the higher 
genus contributions. The singularity of F h (h' < h — 1) does not affect eqs. (|3.5a - 



d). We find that the fixed point condition is the same as eq.( |3.5a| ) except that the 
coefficient bo is now replaced by of the singular term of the genus h free energy. 
The critical exponent condition becomes 

2(1 - h)b h G = b h (2 - r A )(Ai + ^20 ■ 4a 2 ). (3.14) 



The coefficients a% and a 2 of the analytic term are needed to fix these equations. 
However, it is notable that they are precisely the analytic terms in the sphere free 



energy. Therefore the same conditions as the sphere case Q3.6cp and ( |3.6d| ) are 



sufficient to determine them. In this way we find that the fixed point is universal 
for any genus, and the critical exponent Th for genus h is given by 

2-r fc = (l-/i) Tl , 7l + 7o = 2. (3.15) 

This result explains the double scaling behavior for the singular part of the free 
energy 

oo 

F sing (N, g ) = J2 N ~ 2hF L g (g) 

h=0 

= (9 ~ 9*f~ 10 f (N 2hl (g - + less singular terms. (3.16) 
Similarly we observe that the generalization of eq. (|3.15|) to the multi-coupling case, 

2-rf ) = (l-/>) 7l fe) , (3.17) 

holds for each k. 

3.2 Spectrum of continuum theories 

In this subsection we recall the spectrum of operators of continuum theories of quan- 
tum gravity and present the various scaling exponents, to prepare for the following 
subsections. 

Operators of (p, g)-minimal conformal matter coupled to two-dimensional quan- 
tum gravity (hereafter referred to as the (p, q) gravity) derived from matrix models 
can be identified with hamiltonians of Wp-constrained, p-reduced KP hierarchy with 
a source insertion at t p+q |T5], [ll| . Namely the fc-th order KP hamiltonian Ok for 



k < p + q can be identified with the gravitationally dressed (n, m)-th primary oper- 
ator with k = \pm — qn\, 



0\ pm „ qn \ ~ / d l z * n , m e^-^ (3.18) 



\pm— qn\ 

where (p denotes the Liouville field. Operators Ok for k > p + q are the gravitational 
descendants. The W p constraint implies that the free energy 

F(k; U) = jrn h /exp t k Ok) (3.19) 

h=0 ^ fc > X ' h, conn. 

has the following scaling property: 

2 p(k) 

(exp (t k O k + const. O p+q )) K conn ~ t k h {k <p + q), 

(k) . , (fe) (k) 2/c 2{p + q) , v 

r h -To +hj\ \ % = ■ r, 7i = — r- 3.20 

p + q — k p + q — k 

(k) 

The exponent % is related to the more familiar 'gravitational scaling dimension' 



1 A( fc ) by 



A« = l-pk. (3.21) 

2-7o 



These exponents can as well be derived from the Liouville theory; they are related 
with the background charge Q and the Liouville momentum (3 n m of a vertex operator 
by 

^(bm-H) = 2 + _Q_ 7 (bm-gn|) = <3_ 

Pn,m 

/25-c /T. 
3 y pg 




_ V25 — c — \jl — c + 24 Aq 1 '"^ p + g — |pm — gn | 

Pn,m — 



The exponent 7q for the dressed lowest dimensional operator 0i is usually referred 
to as the string susceptibility exponent r str . 



It is instructive to point out that the so-called boundary operators [Pql , Ok for 
fc = mod p or g, should be excluded from the spectrum in our formalism for the 
following reason. Consider the (p, q) gravity as the Wp-constrained, p-reduced KP 
hierarchy with a source insertion at t p+q . Due to the p-reduction condition, KP 
operators Ok for k = mod p decouple 

(O np ■ ■ •) = 0. (3.23) 

On the other hand, KP operators Ok for k = mod q do not decouple automatically. 
In turn they can always be absorbed into the redefinition of the operators of lower 
order than them. For instance, in the case of a unitary theory (q = p+1), shifts in 
the KP time parameters of the form 

ti -> t[ = ti + (polynomial in t 2 , t 3 , • • • , t p _i, t p+1 ) , 
h — > 4 = *2 + (polynomial in t 3 , t p _i, t p+l ) , 



(3.24) 



t' x = t p -i + (polynomial in t p+ i) 



are enough to eliminate t p+ i [^]. In this sense these operators are redundant. This 
type of shift of parameters (of the continuum theory) should automatically be dealt 
with in our procedure to reduce the coupling space, since we have utilized all possible 
reparametrization invariance at the discrete level. Consequently KP operators Ok 
for k = mod p or q would never appear in the spectrum of our RG equation. In the 
subsequent subsections we will observe that this is indeed the case. The fact that p 
and q are treated on equal footing in our formalism should be contrasted with the 
description by the KP flow. This feature will also be discussed in sect. 4. 2. 



3.3 Fixed points of the one-matrix model 



Here we present a full account of the result for one-matrix models, part of which 
has been reported in ref.0]. We investigate the case of two couplings (g^ and g$) as 



an example, by employing the eigenvalue method. The exact solution tells us that 
there exists a pure gravity critical line parameterized by 

2(32^ 2 - z 4 - 6Azw + Az 3 w + QAw 2 ) 16(16z - z 3 - 32w) 

93 ~ z 5 + 2z 4 w - 4z 3 w 2 ' 94 ~ 3(z 5 + 2z 4 w - 4z 3 w 2 ) ' 

= -96z 4 + 3z 6 - 6z 5 w - 6z 4 w 2 + 256z 4 w 3 + 8z 3 w 3 - 128w 4 (3.25) 

and the tricritical point at one end of the line. (See Fig.l). 
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Figure 1: The pure gravity critical line in the one- matrix model. It has the tricritical 
point (0.3066 . . . , 0.02532 . . .) at one end. The dashed curves are unphysical critical lines. 

The dashed curves in Fig.l are unphysical critical lines^|. 

We find four solutions to the set of equations ( |3.10a| — d), as summarized in Table 

1: 



(#3*5 94:*) 


{2- T n 


(p,q), c eff 




(0.3066..., 0.0253...) 


7/3 


(2,5), c eff = 2/5 


O z ~ J 1 e&>w 


7/2 


(432" 1 / 4 , 0) 


5/2 


(2,3), c cff = 


Ox ~ / 1 


-6 


(0, -1/12) 


5/2 


(2,3), c eff = 


Ox ~ / 1 e ft .^ 
Oi ~ / 1 


5/2 


(0, 0) 


-4 


(Gaussian) 




-2 



Table 1: Fixed points of the RG equation for the one-matrix model and 
associated critical exponents. Each fixed point is found to represent a 
continuum (p, q) gravity through identification of operators as listed above. 

* The saddle point equation for the eigenvalue distribution |^j generically has several branches of 
solutions. If we approach the origin within the sheet on which the dashed curve lies, the distribution 
does not converge to the Wigner's semi-circle distribution. Since the identification of the matrix 
model with gravity is based on the perturbative expansion with respect to g around the origin, the 
dashed curves have nothing to do with gravity. 



(2,5) fixed point We find a fixed point (#3*, #4*) = (0.3066 0.0253 .. .) on 
an edge of the pure gravity critical line. Let us describe how this fixed point is 
identified with the (2,5) gravity (Lee- Yang edge singularity coupled to gravity) |f27 . 
By diagonalizing the scaling exponent matrix computed 

_ / 1.5929... -4.5685... \ 
Ui] ~ [ 0.1645... -0.1643... J ' 1 } 

we obtain 2/(2 — ^) = 6/7, 4/7 and Vij = ^ j 4 4722 1" therefore the 

singular part of the free energy is found to behave in the neighborhood of the fixed 
point as 

F sing (iV, g 3 , g 4 )~b- [5g 3 - 4.4722 . . . 5g 4 } 7/3 + b' ■ [5g 3 - 6.2087 . . . 5g 4 } 7/2 + ■■■. (3.27) 

Values of g^ and Gij can be expressed as roots of algebraic equations, though we 
present them in numerical forms for brevity here. 

On the other hand, the content of relevant operators in the (2,5) gravity con- 
sists of Oi, 3 and O5. Among them C 5 is a boundary operator which should 
be suppressed in our formulation. According to eq. (|3.20|) , X ~ / ^ 1)2 e^ 1 ' 2 ^ and 
C 3 ~ / 1 e*' 1 * 3 acquire the scaling exponents 2 — 70 =7/3 and 2 — 7q 3 ' ) = 7/2 
respectively, which are identical to those in eq. ( ^.27] ). This completes the identifi- 
cation, since we generically expect to detect all the relevant operators (excluding 
boundary operators) within the coupling constant space which is wide enough to 
realize the fixed point. 



(2,3) fixed point I We find a fixed point at (g 3 *, g^) = (432 1 / 4 , 0) on the critical 
line. By diagonalizing the scaling exponent matrix computed 

G^( 4 i 5 ^"-2f')/20y (328) 
the singular part of the free energy is found to behave as 



F«ne(N,g 3 ,94) ~ b ■ 



/3I/4 37/4 x 
693 + ( 1 2~) 94 



5/2 



+ b' ■ [^ 4 ]~ 6 + • • • • (3.29) 



This fixed point is identified with the (2,3) gravity (pure gravity) whose sole relevant 
and non-boundary operator is the dressed identity 0\ ~ / 1 e^ 1 ' lV with the suscep- 
tibility exponent 2 — 70 = 5/2. Meanwhile the boundary operator 3 is suppressed 
as expected. 



Here we stress that exponents of irrelevant operators Ou for k > p + q, for which 

(k) (k) 

2 — 7q — 7} < 0, could not be detected correctly within our large- N RG equation. 
It is because identification of iV with the cutoff is based upon the possibility of 



the double scaling limit ( |1.8|) where we let N approach infinity in correlation to a 
combination of ^'s coupled to an operator, sent to a critical value. Obviously this is 
possible only for positive 7^ . Therefore the magnitude of negative exponents such 
as —6 is an artifact of our large- iV RG and untrustable. 



(2,3) fixed point II At (#3*, g^) = (0, —1/12) on the critical line, the G-function 
turns out to develop singularities of the form (3 (gi) ~ O ((VySgj) 5 ^ 6 ), Pi(gi) ~ 



O ((Vij5gj) 1 ^ 6 ) and so forth. These singularities originate from the switchover of 



the dominant solution A to the saddle point equation (|2.30| ) or ( |2.32|) as we pass 



over the g 3 = line |B| . By replacing the G-function in eqs. ( |3.10a| -d) by its regular 
part G reg , we find that (#3*, #4*) = (0, —1/12) is indeed a fixed point with a scaling 
exponent matrix 



/ 4/5 
"I 4/5 



(3.30) 



Since eigenvalues of the scaling exponent matrix are degenerate, the transformation 
matrix V^- in the singular part of the free energy 

FomiN, g 3 , g^-b- [V l3 6g,f 2 + b' ■ [V 2] 5g 3 f 2 + ■■■. (3.31) 

remains indeterminate. This fixed point can also be identified with the (2,3) gravity, 
involving doubling of the dressed identity operator due to the Z 2 -invariance of the 
potential]. 

Gaussian fixed point We find a trivial fixed point at (#3*, g^) = (0,0). The 
scaling exponent matrix indicates the naive scaling behavior 

F sing (N, g 3 , 9i )~b- [5g 3 ]- 4 + b' ■ [5g A ]- 2 + ■■■. (3.32) 



3.4 Fixed points of the two-matrix model 



The two-matrix model Q2.34p has been investigated in ref. |2^|, and is known to 



possess the pure gravity critical line parameterized by 



g 2 = -2cz 2 +zl-c 2 + [c + 4z-y(l + c) 2 + 8cz + 16z 2 ) j , (3.33) 

= (-1 - 2c + 3c 2 + 24c2 + 48z 2 ) 

■(1 + 4c + 6c 2 + 4c 3 + c 4 + 2Acz + 48c 2 z + 24c 3 z + 144cV + 256c,2 3 ) 

and the tricritical point at one end of the line. (See Fig.2). 

^It is well known that the universal part of the free energy derived from matrix models with 
^-symmetric potential is twice as much as that of matrix models without ^-symmetry. 



Figure 2: The critical line in the two-matrix model with the cubic potential. The gravity 
is critical along this line. At one end of the line, the Ising model becomes critical as well. 



We find three solutions to the RG equations in the physical region of the coupling 
constant space {(g, c)\g > 0, c < 0}, as summarized in Table 2: 



(<?*, c*) 


{2- 7 n 


(P,q), Ceflf 


O k ~ f* n , m 


(0.2003..., 0.1589...) 


7/3 
7 


(3,4), Ceff = 1/2 


O x ~ / 1 e^> 1¥3 


(432" 1 / 4 , 0) 


5/2 
< 


(3,2), c cff = 


d ~ / 1 


(0, *) 


-4 



(Gaussian) 





Table 2: Fixed points of the RG equation for the two-matrix model and 
associated critical exponents. 



(3,4) fixed point We find a fixed point on an edge of the pure gravity critical 
line, (g*,c*) = (^(l0(-85 + 62^/19683) 1/2 , -(2y/7 - l)/27^ . By diagonalizing 
the scaling exponent matrix computed 

/ -0.8391... 0.1148... \ 
ij \ -0.0866 . . . 0.3036 . . . ) ' 1 J 

we obtain 2/(2 - 7 J fe) ) = 6/7, 2/7 and V {j = ^ j 'q'^s'" )" Therefore the 
singular part of the free energy is found to behave as | 

F sing (N, g, c) = b- [5^-0.2074 • ■ ■ 5c] 7/3 +b'- [£#+6.3388 ■ ■ ■ 5c} 7 log[5p+6.3388 • • • 5c)+- • 

(3.35) 

* The logarithmic correction should be present in the free energy for the positive integral value 
of 2 — 70 to be detected in the scaling exponent matrix. We have borrowed from the exact solution 
the fact that the correction takes the form 6g 7 (log Sg) 1 . 



This fixed point is identified with the (3,4) gravity (critical Ising model coupled 
to gravity); relevant and non-boundary operators of the (3,4) gravity consist of 
0i ~ Jl e^ 1 - 1 ^, 02 ~ J*i,2 e^ 1 - 2 ^, and 5 ~ J*i i3 e^. According to eq. (|3~20|) , 



they acquire the scaling exponents 2 — 7q 1 ' ) = 7/3, 2 — 7q 2 ' ) = 14/5, and 2 — Jq = 7, 
respectively. Since we have imposed the zero magnetic field condition g + = g_ at 
the discrete level, the dressed spin operator 2 should not be detected. Critical ex- 
ponents corresponding to the rest of the operators (the dressed identity and energy) 
are indeed detected in eq.(|3.35|) and this completes the identification. 



(3,2) fixed point At (g*,c*) = (432~ 1//4 ,0) on the critical line, the c- and a c - 
derivatives of the G-function turn out to develop singularities. By replacing the 
G-function in eqs. (|3.10a| -d) by its regular part G reg as in the (2,3) fixed point II, we 



find this point is indeed a fixed point. Recalling the result of the one-matrix model, 
this point obviously corresponds to the (twice tensor product of) pure gravity. 

Gaussian fixed line We find that all the points on the line g = are fixed points. 
The scaling exponent matrix indicates the naive scaling behavior 

9,c) = b- [5g]- A + b' ■ [5c}° + ■■■. (3.36) 

The existence of the fixed line is a direct consequence of the separation of the 
Gaussian two-matrix model with arbitrary c into two Gaussian one-matrix models. 

Finally it is important to remark that the critical behavior of the second term in 
eq. ([3.27|) , fl3.29|) , or (|3.35|) is realized at each fixed point when one approaches the 



fixed point along the critical line. This means that, by fine-tuning the most relevant 
perturbation by 0\ to zero, we can deform a multi-critical continuum theory while 
keeping the gravity critical. This point will further be discussed in sect. 4. 2. 



4 Renormalization group flow 



4.1 Linear approximation to the RG equation 

In this section we discuss the structure of the RG flow in the coupling space of 
matrix models. For the usual linear RG equation of Callan-Symanzik type 



d 

N-z h7(o) 

dN ny) 



F(N,g)=r(g) + P(g)^, 



(4.1) 



the concept of running coupling constant g(N) is naturally introduced by 

N 



9(N) dg 
go P(g') 



log 



iVn 



(4.2) 



where go = g(N Q ) and the /3-function plays a role of vector field describing RG flow. 



We consider, for instance, the RG flow for the one-matrix model with the cubic 
coupling described by eq. (|2.22|) with g 2 — 1, gz = g, others= 0. Expanding the right 
hand side into a power series in dF/dg around dF/dg\ g= o = 0, we obtain 



F(N,g) 



(4.3) 



Piis) = - 9 2 +3 9 3 J 



9 , o 3 /"T^ff 5 ^ 



[(l-t) 2 + 4^ 2 t 3 (l-(l-^)t)] 1 /2' 



9# 4 , _ g [T^F dt t 



4 



^ 2(1 - g 2 ) + ^ Jo [(l-^) 2 + 4g 2 t 3 (l-(l-g 2 )t)] 3 / 2 ' 

a to) = - 9g5( i: 3g2 + 2g4 W ^ *< 8 



2(1 -# 2 ) 3 * y [(i-t) 2 + 4^ 3 (i-(i-^ 2 )t)] 5 / 2 ' 

Expansion in powers of g shows that faig) ~ —9/2 g 6 + 0(g s ), (3^{g) ~ 9/2 g 9 + 
0(g n ) and so forth. Therefore we expect that the RG flow should be well approxi- 
mated by f3\(g) (~ — g/2 — 3g 3 log g 2 ) for the region \g\ <C 1 where most fixed points 
lie. In fact, we have approximate values of fixed points 

fcfa) = for g* = 0, 0.219345 • ■ • ; fag*) = -1/2, 0.780698 • ■ ■ . (4.4) 

By construction (3\{g) possesses a zero at the Gaussian fixed point with the exact 
value of the slope. On the other hand the deviations of the location of the nontrivial 
fixed point and of the slope from the exact values are 0.237% and 2.472% respec- 
tively. The small magnitude of the errors is obviously due to the suppression factors 
multiplying the nonlinear terms in eq. (|4.3|) as stated above. We emphasize that the 
approximate evaluation of j3 n is useful only after the reparametrization identities are 
taken into account, as is explained in sect. 2.1. 

By the same token, by linearizing the RG equations (|2.33|) around the Gaussian 



values, (dF°/dg 3 , dF° /<9g 4 )| 93=94= o = (0, 1/2), we have also calculated /^-functions 
(G,a 3 (<73j g&] 0, 1/2), G ta4 (g3, g^, 0, 1/2)) for the one-matrix model with cubic and 
quartic couplings. Here we exhibit the real part of the /^-functions in Fig. 3. 



0.1 0.2 0.3 0.4 

g3 

Figure 3: linearized RG flow for the one-matrix model. 
The arrows represent a flow from the UV (N = oo) to the IR (N = 0). 

Similarly the linearization of eq. ( f2l)3| ) around (dF°/dg, dF° /dc)\ g=c= o = (0, 0) gives 
/?i-functions (G tag (g, c; 0, 0), G, ac {g, c; 0, 0)) for the two-matrix model (Fig. 4). 




Figure 4: linearized RG flow for the two-matrix model. 



We immediately observe that they reproduce all the fixed points found in the previ- 
ous section within a very small margin of error. Moreover, IR-repulsive and attrac- 
tive directions of the RG flow around fixed points (which correspond to perturbations 
by relevant and irrelevant operators, respectively) are approximated quite well. 



4.2 Global aspects of the RG flow 



The real parts of /^-functions plotted in Figs. 3 and 4 do not merely reproduce 
the local properties of the fixed points but also qualitatively represent the global 
structure of the RG flow. 

Topology of the RG flow 

We immediately observe from the figures that pure gravity critical lines are 
characterized as the renormalized trajectories emanating from multi-critical fixed 
points. By comparing the exponents with the spectrum of exact solutions (|3.20|) , the 
trajectory of the one-matrix model in Fig.3 is identified with the one corresponding 
to the perturbation of the (2, 5) gravity by its dressed (identity) operator, and 

that of the two-matrix model in Fig.4 to the perturbation of the (3, 4) gravity by 
its dressed ^(1,3) (energy) operator. Both of them flow into the pure (2, 3) gravity. 
These results lead us to propose the following conjectures: 



perturbation of an UV theories of two-dimensional gravity by its least relevant 
operators leads to the IR theory with the neighboring order of criticality; 

the effective central charge^] c^ q ^ = c^ p,q " > — 24A™ in = 1 — 6/pq decreases along 
the RG flow from the UV to the IR. This c e fj counts the density of states of 



the matter sector [Eg 



In particular, the RG trajectory corresponding to the dressed ^(1,3) perturbation 
interpolates two neighboring unitary models. This is an intriguing example of the 
gravitational counterpart of the familiar result on two-dimensional renormalizable 
field theories [13j . One may naively anticipate this result from the knowledge about 



the statistical and field theories over fixed backgrounds. However, we have vanishing 
total central charge of the system when we couple the matter with quantum gravity. 
Therefore it is nontrivial to establish the concept of 'reduction of degrees of freedom' 
along the flow from the UV to the IR in presence of gravity. 

We remark on the RG trajectory between two (2,3) fixed points in Fig.3. At first 
sight it is curious that a RG flow apparently exists to interpolate two distinct fixed 
points (the UV (2,3) fixed point II and the IR (2,3) fixed point I) corresponding to 
the same continuum theory. We interpret this fact as follows. The one-parameter 
subspace (#3 = 0, #4) is stable under the RG flow, because our RG transformation 
maintains the Z2-invariance <fi — > — 0. Moreover, the exact solution shows that the 
singular part of the free energy is discontinuous as we trace it on the critical line from 
the region g 3 > to g 3 = 0. As the g 3 -derivatives of the free energy is ill-defined, we 
can attach no meaning to an RG trajectory connecting the Z 2 -invariant subspace 
(?3 = with other regions; we can make sense only for the eigenvector of the scaling 
exponent matrix and the RG trajectory which are along the g± axis around the (2,3) 

§ The effective centrai charge in a conformal field theory is defined as the measure of the rate 
of divergence of the canonical partition function of the theory 

Z((3) = Tr e _/Hr ~ f3 r exp(7rc off /6/5) (H = 2tt(L + L - c/12)) 

as (3 — ► 0. Since Z((3 — > 0) is interpreted as the regularized form of the number of states in the 
theory, c c s describes the asymptotic growth of density of states. 



fixed point II. 



Comparison to the KP flow 

Let us further contrast these results with the integrable hierarchical description 
of two-dimensional gravity We again recall the fact that (p,q) gravities are 



universally described by the Wp-constrained, p-reduced KP hierarchy with a source 
insertion at t p+q (U denotes the time parameter coupled to the z-th order Hamiltonian 
Oi of the hierarchy). We list the sets of relevant operators (usually referred to 
as gravitational primaries) in (3, 4) and (3, 2) theories in KdV description of the 
matrix models as well as the corresponding operators in the Liouville field theory 
description: 

(3,4) KP operators: O x 2 ■ (0 4 ) C 5 • [O r ] 

DDK operators: J* (1>1) /* (1>2) f%, 3) 

i 

(3,2) KP operators: O x (0 2 ) ■ (C 4 ) [O s ] 
DDK operators: J %i) 

Table 3: relevant operators of (3, 4) and (3, 2) theories 



In this table supplemented are two types of KP operators: 3+q corresponding to 
the critical action S^ 3 '^ parenthesized by [ ], and boundary operators Ok, k = mod 
q parenthesized by ( ). It is naturally expected in the KP description that the (3,4) 
gravity perturbed by 5 should flow into the (3, 2) gravity; since 5 is relevant, t 5 
will increase to infinity along the RG trajectory against a fixed source value of t 7 . 
Then by rescaling each by {t^)~ 1 ^ so as to render £5 finite and t 7 infinitesimal 
while keeping the form of the hierarchy, we can describe the IR theory by the same 
3-reduced KP hierarchy having a source insertion at £5. This theory is regarded as 
the (3, 2) gravity. Our result on the two-matrix model is in accord with this picture. 

On the other hand, the sets of relevant operators of (2, 5) and (2, 3) theories are: 

(2,5) KP operators: O x ■ O a ■ (C 5 ) • [0 7 ] 
DDK operators: /*(i, 2 ) ^ (2 ' 5) 

\ 

(2,3) KP operators: O x ■ (C 3 ) • [O s ] 
DDK operators: / 

Table 4: relevant operators in (2, 5) and (2, 3) theories 



It is awkward to interpret the RG flow within the framework of the KP hierarchy 
alone, as we describe below. One naively expects that the O^ operator perturbation 
would give the (2, 3) gravity from the (2, 5). However, the boundary operators such 
as O5 is physically redundant, since it can be expressed by means of other primary 
operators. In fact we have found that the RG flow from the Os (~ / 1 e^ 1 ' 1 ^) operator 
perturbation of the (2,5) gravity leads to the (2,3). Our large-iV RG treatment, 



which suppresses appearance of such redundant operators automatically by the use 
of reparametrization identities, is efficient indeed in such cases for identifying the 
RG trajectories and perturbing operators. 

We remind the reader that the 'universal' description of (p, *) theories by p- 
reduced KP hierarchy is derived by blowing-up one of the critical regions in the 
bare coupling space {gj}- In other words the KP time parameters {ij} coupled to 
continuum operators are renormalized coupling constants defined with reference to 
a specific critical point. Hence it is a priori not guaranteed that the p-reduced KP 
hierarchy can describe an RG flow interpolating between two distant critical points 
(p, q) and (p, q') by a limiting procedure with respect to £j. 



5 Discussions 



In this paper we have constructed the RG equation for matrix models by regarding 
N as the cutoff, and investigated its implication. First we have performed the 
RG transformation N + 1 —>■ N to obtain the RG flow in the enlarged coupling 
space. We have reduced it into the original finite dimensional coupling space by the 
use of reparametrization identities, at the cost of allowing the nonlinearity of the 
RG equation. These reparametrization identities are the discrete Schwinger-Dyson 
equations for the sphere correlators (|2.9| ) and ( |2.49|) . It is crucial for this procedure 
that these equations close algebraically to determine the required resolvent. The 
RG equation is found to give the critical points of the exact solution. It also exactly 
reproduces the spectrum of relevant operators of the Liouville gravity at its fixed 
points. After confirming the adequacy of linear approximation, we have investigated 
the RG flow in the coupling space. From these examples, we conjecture that a kind 
of the gravitational analogue of Zamolodchikov's c-theorem holds in accordance with 
some previous expectations [17 . 

One of the characteristics of our treatment is that it suppresses the boundary 
operators which inevitably arise from the usual KP description of quantum gravity 
derived from matrix models, but are absent in the BRST analysis of Lian and 
Zuckerman's [p9f] . A direct consequence of this is that we do not need them to 
unambiguously determine the RG flow. Further, if applied to multi-matrix models 
it can in principle describe the RG flow which changes both p and q of the (p, q) 
theories that is impossible to realize in the KP description. Work in this direction 
is in progress. 

In order to establish the 'gravitational c-theorem', it is necessary to construct a 
counterpart of the c-function in the framework of matrix models. The c-function is 
a potential function for the /3-function, which should monotonically decrease along 
the RG trajectory from the UV extremum to the IR one. It might exist, since the 
G-function itself plays the role of the c-functions in nonlinear RG equations. 

As originally motivated Brezin and Zinn- Justin, another intriguing application 
of our method is to solve matrix models which do not allow angular integration. 
Such models include candidates for a system of c > 1 matter coupled to gravity. For 
instance, the n-Ising model over randomly triangulated surfaces can be described 



by a matrix model with 2 n -plets of N x N matrices 0i, • • • , 02™ as: 



Z(N,g,P) = JU^^r expi-NtiVfa, •••,</>*>)} 



V((f)i, (f> 2 n) = -^0 r A rr /0 r / + gZl^r 

t,t' t 

t = (01, cr 2 , • • • , a n ), Oi = ±1, 



V fa) 

It is straightforward to write down the RG and saddle point equations, 



NjL + 2 
dN 

I 1 



F(iV,^) = (itry^..,^) + -i 



(as) 1 



v9 



+ ^Tr log (1 A rT ' + g ((a s ) 1<W + </>AtO) ) - 2"" 1 + O y y J . 



P 2 
1 



(5.1) 



+ 2 n | (a s ) 3 



(5.2) 



1 



(2 cosh/?)" 



(a B ) +g(a s ) 2 + g 



N 



1 



,lA + 5 ((a B )l ( J + 0) 



(5.3) 



where the trace over the (2 n ■ AT) x (2™ • AT) matrix is denoted by Tr, and the trace 
over N x N denoted by tr. We have assumed a SjT does not depend on r in the same 
way as in the analysis of the two-matrix model. Thus the problem reduces to the 
computation of the resolvent 



W(z) 



1 



:Tr 



1 



N zl5 + (lA + g<f>) 



(5.4) 



by using the reparametrization identities. So far we have not yet seen whether the 
resolvent can be determined by a closed set of Schwinger- Dyson equations or not. 
Even if we do not obtain a finitely closed set, however, we can exploit available 
approximation schemes, such as the linearization of the RG equation. 

(Non-polynomial) nonlinearity is inevitable in constructing an RG equation obeyed 
by the free energy F(N,g). Consequently nonlinearity makes the identification of 
fixed points and exponents complicated, and makes the concept of running coupling 
constant unclear. While we are writing the manuscript we noticed two recent papers 



which may be relevant to this point. In ref.|30[ the free energy of the Penner matrix 
model is shown to obey a linear RG equation. It is of great interest to examine 
whether the linearity of the RG equation is the peculiarity of the Penner model 
or a generic feature in a certain class of matrix models. A more recent paper J3T 



suggests an interesting possibility for constructing a linear RG equation by using a 
mapping of matrix models to 0(A^ 2 )-vector models. It was also proposed to con- 
sider the RG equation for the first derivative of the free energy dF/dg. Although 
generalizations of this mapping to the subleading orders in A^ are hard to provide, 
so far we have observed that a certain simplification occurs when we perform an RG 



transformation for dF/dg; there appears neither the tr log-term present in eq. 



nor the integration term in eq.( 2.22|) . 

We hope to report on these topics in a future publication. 
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A Loop equation for the two-matrix model 

In this appendix we exhibit the explicit form of the one-point function of the resol- 
vent operator (Wq(z)\ in eq.( |2.50| ) by using the reparametrization invariance of the 
partition function ( |2.34| ). We will consider the following reparametrizations: 

cf>' = </> + e Pa (za + -a x + <f>\ P a , (A.la) 

$ = </> + e £ p a jo-!, (^ Z a + -o"! + ^ | P a , (A.lb) 
0' = <j) (A.lc) 
+ e P * (^ ztr o + + <t^j + [zct q + -^T X + cr^j P a , 

0' = (f> + tY,Pa jcn</>cri, ^cr + -^j x + cj^j I P a , (A. Id) 

These reparametrizations maintain <j)± hermitian. Invariance under these reparametriza- 
tions induces the following reparametrization identities 



\(w^)) 2 +{gz 2 -z+ c - g 
1 1 







) (WoM) - 







+c (W n (z)) + g (^Tr<^ + 2 - 2gz, (A.2a) 
= \(Mz)){Wi{z))+(cz-- g - C ^\(^ 



-c(Wu(z)), (A.2b) 

= -^(# w) 2 -^(#oW)(^ 1 (z)) + ^-c^(iy (z)) 

+ (-1 + z 2 - gz 3 ) (Wtiz)) + f-~ + cz^j (W n (z)) - c (W wl (z)) (A.2c) 



where 



+c(W 111 (z) 



\ (W n (z)) (W (z)\ + - gz) (w (z)\ + -z(l - gz) (w x {z] 
£ 9 9 



-z(l - gz) {W n (z)) + c(W m (z)) -c(W ni (z)) + (1 - gz) ( -Tr0 



+g ( — Tr<ri0<ri0 



w hh-h( z ) 



— Tta j x (j>a j 2 (j> ■ ■ ■ <pcr jk ( zer + -(T 1 + 4> 



N ' 



9 



(A.2d) 
(A.3) 



After some computations, we obtain a quartic equation obeyed by {W (z 



z + cz + gz A 

9 



\ (W (z) 

~(1 + 2c - gz) (W (z)) 2 + (2.% + ^--z)(c + gz)(l-gz) (W (z 



2g 

+ c{l - gz) l^—Ticf) 



g 2 g 2 g 



2gc(—Ti(T 1 (f)cr 1 (p 



\(w (z)) 2 + 



~ CZ ~\2g + Z ) (1 " 9Z) 



W (z))-c 



(A.4) 



+ -(l + 2c-gz)(c + gz)(l- gz) 

g 



/ 1 



+ cz + gz 2 1 (W (z)) + g (^Tr0) +2-2gz 



There appear one-point functions (tr0 + 0_) = (1/2) (Tr<7i0<Xi</>) and (tr(0 + + 
(/)-)) = (Tr0) in the above expression. We express the latter in terms of (tr0 + 0_) 



and ( tr (<jf + + 



by solving the first two equations in ( |2.49|) , 




1 



I 1 



x ^tr0 ± ) + ^(-tr^)+c(-tr0 T 



/ 1 



\N 



N 
I 1 



/ 1 



-tr«/4)+^-tr<^)+c(-tr<M>_). 



\N 
/ 1 



\N 



(A.5a) 
(A.5b) 



Promoting dF/dg = ((l/3iV) tr ((f> 3 + + <p 3 _)) and dF/dc = ((1/N) tr to inde- 

pendent variables, eq.( |A.4j ) with this replacement determines the one-point function 
of the resolvent 

/ OF dF\ 

(W (z)) = W lz;g,c;—, — \. (A.6) 
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